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,^ , Abstract: Using recently proposed method of discrete Hirota dynamics for integrable 

(1+1)D quantum field theories on a finite space circle of length L we derive and test 
^ numerically a finite system of nonlinear integral equations for the exact spectrum of energies 

f~^ ■ of SU{N) X SU{N) principal chiral field model as functions of rriL, where m is the mass 

r^ ' scale. We propose a determinant solution of the underlying Y-system, or Hirota equation, 

in terms of determinants (Wronskians) of N x N matrices parameterized by A^ — 1 functions 
r — ^ . of the spectral parameter 9 with the known analytical properties at finite L. Although the 

^P i method works in principle for any state, the explicit equations are written for states in the 

C/(l) sector only. For A > 2, we encounter and clarify a few subtleties in these equations 
related to the presence of bound states, absent in the previously considered A = 2 case. 
L^ In particular, we solve these equations numerically for a few low-lying states for A = 3 in 

;_, ■ a wide range of rriL. 
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1. Introduction 

Integrable 1 + 1 dimensional quantum field theories on a finite space circle were rather 
intensively studied in the last 20 years |||, ^, ^ ^ |5|, |^ |^, |^, ^, ^] . A great deal of success 
in the exact treatment of the finite size effects in various integrable QFT's was due to the 
thermodynamic Bethe ansatz (TBA) approach resulting in a system (in most of the cases 
infinite) of non-linear integral equations. It was realized that the TBA equations could be 
rewritten in a functional, Y-system form 0. 

Recently, a novel, and hopefully quite general approach to this problem was proposed 
in |1^ based on the integrability of the Y-system. The Y-system is known to be a gauge 
invariant version of the famous Hirota bilinear equation, often called the T-system, in 



its discrete form|12]. The underlying discrete Hirota dynamics is integrable and general 
solutions with various boundary conditions, appropriate for various problems, from matrix 
models to quantum spin chains and quantum sigma-models, can be found in a determinant 



(Wronskian) form [13|. For finite rank symmetry groups the matrix of the Wronskian 
contains only a finite number of functions of the spectral parameter, drastically simplifying 
the problem: the infinite Y-system is reduced to a finite number of non-linear functional 
equations for these functions. Then the subtlest point comes: We should guess the analytic 
properties of these functions w.r.t. the spectral parameter. This is relatively easy to do 
for the spin chains where the polynomiality of transfer matrix leads to the final answer 
in terms of a set of Bethe ansatz equations. For the QFT's at a finite volume L (length 
of the space circle) the situation is much more complicated and the analyticity properties 
of the Y-functions are not so obvious. Nevertheless, it often appears to be possible to 
extract them, partially from physical considerations, partially from certain assumptions of 
the absence of unphysical singularities. It helps to transform the Y-system into a system of 
non- linear integral equations (NLIE), more tractable, and better suitable for the numerical 
studies. The resulting equations can remind the Destri-DeVega NLIE P] or even coincide 
with them for the (rare) models where these DDV equations are known. 



This program was performed in |11] for the SU{2) x SU{2) principal chiral field (PCF) 
for a general quantum state, and the numerical study of the finite size spectrum was 
successfully done for a variety of interesting states, from the vacuum and mass-gap to quite 
general states, in the so called U{1) sector or even lying out of it (i.e. having excitations 
in left and right SU{2) modes). 

In this paper, we will construct within these lines the corresponding NLIE's for SU{N) x 
SU{N) PCF at any N. We use the Wronskian solution of ||l^ for the underlying Hirota 
equation in terms of determinants of N x N matrices and guess the correct analytic form of 
the functions entering the Wronskian. For the vacuum state, the asymptotic Bethe ansatz 
(ABA) based on the scattering theory and, strictly speaking, valid only for sufficiently large 
length L, teaches us that there are no singularities on the physical strip of the rapidity 
plane, at least for not too small L's.^ 



^The argument can't exclude the possibility that at sufficiently small size, some extra singularities occur. 
However, our numerics give serious evidence that at least at A'^ = 3 such extra singularities do not appear. 



For excited states there are certain poles entering the physical strip, and their qualita- 
tive structure is guessed from the ABA. The explicit construction is done only for states in 
the U{1) sector, but we sketch out the generalization to any state. We show how the exact 
S-matrix of the model (including the CDD factor) naturally emerges from this approach 
by simple analyticity assumptions. 

The presence of additional singularities on the physical strip related to the bound 
states, absent for N = 2, leads in A^ > 2 case to significant modifications, already in 
the expression for the energies of excited states. We find from our NLIE's the finite size 
(Liischer) corrections which reveal the presence of the so called //-terms. We also test our 
NLIE's analytically, comparing the results with the known analytic data in the ultraviolet 
(conformal) limit. Finally, we demonstrate the power of our approach by solving the 
resulting NLIE's numerically, for the vacuum energy and the energies of some low lying 
excited states as functions of the size mL for A^ = 3. 

One of the principal motivations for our work is the hope to fulfill the same program in 
the case of recently constructed AdS/CFT Y-system [|l4| , [l5| , |lq ] for the exact spectrum of 
anomalous dimensions in N=4 super symmetric Yang-Mills theory. The PCF model, having 
A — 1 particles (including A — 2 bound states) in its asymptotic spectrum, bears many 
similarities with the AdS/CFT case where the number of bound states is infinite. The 
corresponding Wronskian solutions of AdS/CFT Y-system, or Hirota equation with the so 
called T-hook boundary conditions is already available |17, O]. 



2. The principal chiral field model in the large volume 

In this section we will give the definition of the PCF model, remind the reader the scat- 
tering theory for the physical particles and the ABA equations, and describe the finite size 
equations in terms of the Y-system. 

2.1 The PCF model, its S-matrix and the large L ABA 

The SU{N) X SU{N) PCF model has the classical action 

5pcF = — \ f dt f dx tv{h-^dahf , heSU{N). (2.1) 

2eo J Jo 

The spectrum of this asymptotically free theory in the infinite volume L — )■ oo consists 
of A — 1 physical particles with masses 

sin ^ , , 

ma = m—^ 2.2 

sm^ 

where the lowest mass scales with the bare charge eo according to the asymptotic freedom 

47r 

m = —e ^"o ( A is a cut-off). Its wave function transforms in the fundamental representa- 
tion under each of the SU{N) subgroups. The exact S-matrix for bi-fundamental particles, 
found from the conditions of factorizability, crossing, unitarity, analyticity and the bound 



state structure [|l9| , reads^ [§0|, |2l| : 



SuiO) = XcnAO) ■ 5o(^)^^ Soie)^^ 



(2.3) 



SoiO) 



r(-4)r(i^) 



Xc 



smh(7r6'/iV + in/N) 
smh{TTe/N - in/N) 



(2.4) 



where we introduced the standard SU{N) R-matrix RL,R{d) = O + zPl^ji and P is the per- 
mutation operator exchanging the left/right spins of the scattering particles. In particular, 
crossing and unitarity lead to the following identity 



JV-l 
2 



;Ar-l 



\N' 



n So{9 + zk) = {-iy ^^^^ 



(2.5) 



JV-l 
2 



on the scalar (dressing) factor. 

We can use this S-matrix to study the spectrum of A^ particles on a periodic space 
circle of a sufficiently big circumference L ^ m~^ imposing the wave function periodicity 



N k~l 

j=k+l j=l 



(2.6) 



which quantizes the momenta of the physical particles. The asymptotic spectrum is then 
given by 

E^^mcoshi—6j]+0{e-"'^) (2.7) 



where 9j are solutions to the system of nested Bethe equations following from the diago- 



nalization of (2.6). This diagonalization can be performed by means of the algebraic Bethe 
ansatz and leads to the ABA equations ( 4.20| ) and (4.9) 



, 23|.'^ We will derive these 
equations as a large L limit of the Y-system of the model - a system of equations valid at 
any finite volume L and presented in the next subsection. 

The eq. ( 4.20| ) represents the diagonalized version of the periodicity condition (^.6|). 
(|4.9| ) is the set of 2(A^ — 1) nested Bethe equations for the auxiliary right and left magnon 



(k) (k) 

roots u, and v, following from a regularity condition, as we will see in section |4J. 



Note that the ABA equations ( [4. 9] ) remind the Bethe ansatz equations for two inhomo- 
geneous SU{N) spin chains with the inhomogeneity parameters 9j given by the rapidities 
of physical particles. Their dynamics is defined by the periodicity equation ( 4.2[1| ). So the 
large L limit can be also called the "spin chain limit". 



^In the N — 2 case, these definitions give Xcdd ~ ~li which can also be expressed by multiplying So 
by i. 

^In what follows we will measure all dimensional quantities in the units of the mass m, so that we put 
everywhere m — 1. The only continuous parameter of the problem is now the volume L. 



I" 




0.0 



Figure 1: The (a, s)-strip for Y-system and T-system 



2.2 TBA, Y-system and Hirota equation 

The generalization of these ABA equations to any length L is achieved by the TBA trick 
|2|: the system is put on the space time torus, with a finite space period L and a big 
euclidean "time" period R ^ oo. Then, using the relativistic invariance, we exchange the 
roles of time and space and rather solve the problem for the same system but with infinite 



space extent R placed into a periodic time L interpreted as the inverse temperature |24|. 
The full energy spectrum of such an infinite system can be found from the nested BAE 
(|2.6|) and from (2.7) by means of the so called string hypothesis. The resulting equations 



for the densities of bound states are presented in ||25|] , following the direct solution of the 
PCF given in p^ (see also p3l)- "^^^ free energy calculation at a finite temperature for 
such an infinite volume system can be done thermodynamically, using the saddle point 
approximation due to |£^]- Then the resulting integral TBA equations can be rearranged 
into the Y-system'* 

1 + y. ,+1 1 + Ya,.-1 



-" a,s -' a,, 



1 + (Ya+1,.)-1 1 + {Ya-l,s)-' 

a = 1,2, . . . , N — 1; — oo < s < cx3 



(2.8) 



where, by definition, lo,s 



Yfyf s = oo and we have the following boundary conditions at 



J^a,, 



_ g-Lp„(0)5,.o X constat. 



Pa 



, ^26l7r sin( 
cosh(— — 



OTT ^ 

N J 



N 



siHn) 



(2.9) 



This Y-system describing PCF is an infinite set of functional equations (2.8) with the 
functions Ya^si^) of the spectral parameter 9 defined in the nodes marked by black and 
white bullets in the interior of the infinite strip in a, s lattice represented in fig. ^/\\ . 

A direct but rather tedious derivation of this Y-system was performed in the Appendix 
A of |11| for N = 2. The generalization of this calculation to A^ > 2 is rather straightfor- 
ward, but the Y-system ( p.8D is known from other considerations [29] and is a very universal 
system of equations describing the integrable Hirota dynamics. 

*To make many formulas less bulky, the shifts of the spectral parameter will be often denoted as follows 
/± = f{9 ± i) , /±± = f{e ± i), and in general /'±'=1 = /(^ ± ifc). 



As we will see later, the expression for the momentum Pa{0) is the only one compatible 
with the y-system and relativistic invariance, up to a normalization that can be absorbed 
into the definition of the size L of the spin chain^. As a result of ( |2.9| ) we see that the 
middle node Y-functions, Yafi, a = 1, 2, • • • , A^ — 1, are exponentially suppressed at large 
L or at large \0\. 

Obviously, the Y-system (2.8) has many solutions and to specify a solution we have 



to describe its analytic properties. To have a qualitative idea of the analyticity we have 
to consider a certain limit for the solution where we know the corresponding Y-functions 
entirely as analytic functions of the spectral parameter 9. The most convenient limit is 
L — )• oo where we can solve the Y-system directly, with the appropriate physically natural 
analyticity assumptions, to obtain explicitly all Y-functions and make a link with the exact 
scattering matrix and the resulting ABA equations, as it was done for the N = 2 case in 



1 11]. We will give this asymptotic solution in the next section. Then the Y-system, in the 
form of TBA equations, can be in principle solved numerically by iterations, starting at 
large L and then adiabatically approaching finite, and even very small L's corresponding 
to the ultraviolet CFT behavior. The method was successfully used for various integrable 
sigma models, including the SU{2) PCF |11, 30, ^, |^. It will be also the main method of 



this paper devoted to the SU{N) PCF for N > 2. For the vacuum state, the information 
from ABA is trivial: it suggests that we don't have any singularities in the physical strip 
—iN/A < Im(^) < iN/4:, at least for not too small L's, since there are no Bethe roots. ^ 

The TBA procedure described above leads to the following expression for the vacuum 
energy 

N 



m ^ sinf — 1 r°° /97r \ 

E^acuumiL) = -^ V -4fY / d^cosh —9 log (1 + Ya,o{9)) (2.10) 



With certain modifications in the analytic properties of Y-functions, described in the 
next section, the equations (2.8-2.10) appear to be appropriate not only for the vacuum 



state, as it was originally derived from the string hypothesis, but also for the excited states 
P, 0]. y-functions for various excited states differ by their analytic properties which can 
be qualitatively inferred, as it was mentioned above, from the same states in the ABA. 
A naive heuristic proposal which worked well for N = 2 case is that the excited states 



correspond to the appearance of logarithmic poles in the integrand of ( 2.10 ) at the points 
9j where 

Yi^o{9j + iN/4) + 1 = (2.11) 

If the contour is deformed so that it encircles these singularities, the pole calculation will 
give a contribution X^,=i w-cosh ^^ which fits well the prediction of the ABA formula 



( |2.7| ). However, the situation appears to be more complicated at A^ > 3, already because 
of the fact that unlike the N = 2 case of [^], the solutions 9j of ( p. 11 ) are not necessarily 



real and this naive prescription should be slightly modified in order to get a real energy. 



^So that the length L is actually measured in units of mass. 

®It does not guarantee that we will not have some singularities entering the physical strip when L 
becomes small enough. But our numerical result don't suggest such a strange behavior. 



This will be explained in detail in the section |5|. One should admit that the right formulas 
for energies of excited states in the integrable sigma-models are still rather a matter of a 
natural guess then of a reliable derivation. More insight is needed into this issue. 

To solve the y-system equation (2.8) we will often use it in the form of the Hirota 
equation 



^a,s^ a,! 



Ta+l,sTa-l,s + Ta^s+lTa,s-l 



(2.12) 



on a set {Ta^s} of functions of the spectral parameter 9 related to the original y-functions 
as follows 



Yr 



a,s 



The Hirota equation ( p. 12 ) is invariant under the gauge transformation 

rp , [a+s] [a-s] f-a+s] f-a-slm 

Tns -^ x\ Xo X^ Xa Tas 



(2.14) 



so that T-functions are gauge dependent, whether as y-functions (2.13) are gauge invariant. 
Another useful relation following from ( |2.12 ) is 



l+Ya. 






(2.15) 



3. Central node equations 

The central node Y-functions Yafi related to the black, momentum carrying nodes on the 



fig. |2.1| play a special role in the Y-system. It will be useful for the future to solve the 
corresponding Y-system equations for these functions entering the l.h.s. of (|2.8| ) in terms 
of the r.h.s. 

Let us rewrite the y-system (E]q) in the form 



Y+Y~ 



1 + Ya.+l 



1 + Y.^ 



a-s„ 



^l-<5a 



\l-<5„ 



(y„+i,,)^-''-~- (ya-i,.)^-"-^ (1 + Ya+i,sy 

At s = it can be rewritten using ( |2.15 ) as follows: 



[l + Y, 



a-s„. 



(3.1) 



a— 1, 



y 



*A 



a,0 



rp*A rp*A 

-'a+l,0-'a-l,0 



rp+ rp- 
rp rp(L) 



rp rp(L) 

.-'0,1^0,-1- 



(3.2) 



where we introduced a discrete D'Alembert operator A on the interval a G [1, A^ — 1] defined 
by the formula 



K 



i*A 



F+F~ 



^l-Sa 



\l-<5a, 



(3.3) 

for any function Fa{0), and where 6 is the Kronecker symbol, used to add the counter terms 
necessary to satisfy ( |3.1| ) even at a = 1 and o = A^ — 1. By a subscript (L) in T^ J_-^ in 



''The terming "discrete D'Alembert operator" becomes clear if one takes the logarithm of the r.h.s. and 



the l.h.s. of (3.3) 



we denoted a T-function in a gauge of a type {L) which can be different from the gauge of 
the other T- functions in that formula. We can do that because 1 + Y„ ^ = , °'x in 

-'O.s ^N.s 



the right hand side of (3.1) are gauge invariant, and we are allowed to write each y-function 
in terms of T's each taken in a different gauge. The meaning and the notation of the gauge 
(L) will be explained later. 

We can act by A~^on both sides of (|3.2D, to get 



YaA 



'Lpa{e)_ 






rp+ rp- 



rp rp{L) 

.^Cl-fo,-!, 



<Sa,l 



*A- 



(3.4) 



where pa = cosh(-^) sin(^). The first factor is the only zero mode of A fitting the 
relativistic invariance and it gives the right large L asymptotics (p^). Furthermore, the 
action of the operator A~^ can be easily identified by the fourier transform in 6, a variables, 
so that the final expression is ^ 



Yr 



a,Q 



-Lpa(e)^ 



rp rp{L) 



T, 



T„-i 



n 



N-a 



a+lfl-La-lfi 



rp+ rp- 
rp rp{L) 



n. 



rp+ rp- 

rp rp{L) 



*K^ 



(3.5) 



where the "fusion" operator n^ is defined as the following product 



{k-l)/2 

n /(^+^i) 

j=-(k-l)/2 



and the kernel Kpf is the operator inverse to Hn- V/ regular, [Ui^lf]) 
transform is 



*Ka 



1 



E 



JV-l 
2 



(3.6) 

/. Its fourier 

(3.7) 



iV-l 
2 



p2injuj 



Back in the 0-space it takes the form 
Kn 



tan (ivr (^ 



2ie 

N 



))+tan(l.(^ + M)) 



2N 



(3.8) 



4. The large L, "spin chain" limit of Y-system and its relation to ABA 

We will derive in this section the large L, ABA equations ( [4.2Cl| ),(^^) directly from the 
Y-system (^^). Following the logic of [^ we use the fact that the Y- functions of the 
momentum carrying (black) nodes are exponentially small in this limit: 



y 



This implies that the two wings, left (for s < 0) and right (for s > 0), of the Y-system 
p!q ) are almost decoupled and can be treated separately. 



*We use the notation /*^" for e'"^-**^" where * stands for convolutfon. 



4.1 Expressions for T-functions in the spin chain limit 

Eq.(^) suggests that either T^^i ~ e'^P^C') or r^^i ~ e^^^^^^l Which one does so 
(whereas another one is finite) is a matter of choice of a gauge for T-functions. 

We will work with two different gauges (R) and (L), such that in the large L limit we 
have 

TJ.% « 1 , ^i? « 1 , T^%o - 1 > ^i?<o - 1 (4-2) 



l<a<N-l 

(R) will be called the "right-wing-gauge" and (L) the "left-wing-gauge", and when this 
superscript will be omitted it will be implicitly assumed that we are working in the (R) 
gauge. 

In the large L limit, the T-functions of the left (L) and right (R) gauge both describe 
the same Y functions but (up to exponential corrections) they satisfy Hirota equation 
restricted to the wings s > (resp s < 0). Moreover, these T-functions are in this limit 
analytic on the whole complex plane, and therefore polynomial. 

Such a solution of Hirota equation is well known in applications to the fusion pro- 
cedure in similar spin chain systems, bosonic [O] or even super symmetric |33, R3]. First 



we parameterize Ti^g in terms of N functions XK (9), j = 1, . . . , N hy means of the 



following generating functional 



^(j) 



-1 / Mjn .^ \ -1 



VrW = (l-X(%)(^)e*^«)" {l-Xl^l^^{e)e^^<>y ...{l-Xl^\9)e'^^y (4.3) 

(4.4) 



- r(J)(^ + §(.-!)) 



These functions can be further expressed as follows 



X(^) - ^^ ^ A- - 1 2 TV r4 5^ 

in terms of some Q'^'-functions, where the superscript W = R,L indicates the wing 
that we study (either right or left). These Q functions in the corresponding gauge are 
polynomials characterizing different solutions of Hirota equation in the large L limit, their 
roots are the Bethe roots describing various excited states - solutions of the y-system^: 





- .f 


AT 


qH 



Qf^'-He)^! (4.7) 



^We assume here that there exists a gauge such that the large L hmit is described polynomial Q[, 
functions. Although it needs a better understanding from the point of view of Y-system, it is the case if we 
start treating the large L limit from the S-matrix by the ABA approach. 



In particular, we have from p.3| 



N 



TrW = V^(^-f)E4TW (4-8) 



fc=i 



Ti^i should be free of poles, i.e. polynomial. But for each Bethe root Wj = u- or Wj = v • , 
the two functions X.^,-, and X.^J^. have a pole at the same position Wj — § ("2" ~ ^) • -^y 



requiring their cancellation in the sum (|4.q ), we get a constraint on the position of Wj, 
which we will call the auxiliary Bethe equation: 



, _ QJT K- - V2) Qf /^V, + ^)Qfl^\w, - ^/2) J A: = 1, 2, . . . , iV - 1 

gf_/^) (^, + V2) Qf /^^ {n,, - ^)Qi%'^ (u., + ^/2) ' j -, = -f ^ resp vf^ ^ ' ^ 

The rest of the T-functions in the right wing can be expressed through the Bazhanov- 



Reshetikhin determinant [34| 



_ deti<j- fc<a Ti^s+k-j {e+^{a + l-k-j)) 



and they are also automatically polynomial in virtue of (14. 9| ). 

Among these Q-functions, the polynomial function Qjy = cp, encoding, as its roots, the 
rapidities of all physical particles, will be of a particular importance, and the vanishing of 
7a -1 o^ ^a 1 "i^^ ^o (|4.1|) implies the following asymptotics^*^ 

^(^)= limrij^)(e + ^^^^) (4.11) 

= lin. T<Ji„(« + .^) = Inn t'^ „(« + ,^) (4.12) 

L— s>oo 4 L— s>oo ' 4 

= m. ri^i„((- - i^) = lim li^.>,„(« - i^) (4.13) 

These relations translate all the zeroes 9j of Yifi{9 + i^) + 1, according to the Bethe 



equation ( 2.11 ), into the zeroes of T-functions. 

4.2 Asymptotic Bethe ansatz (ABA) 

Now we will reproduce from the Y-system in large volume limit L ^ oo the ABA equations 



(4.20),( |4l^ ) for the spectrum of energies. In this spin chain limit, ( p.5| ) can be employed to 



compute Yafi to the leading order, by using the asymptotic behaviors (4.2). 



At this point, it is interesting to notice that the crossing relation implies that up to a 
zero mode of Kj\[ 

(p+\*^^ J+2-N] / [-2N+1]Y^^ [+2-N]g[~N] [+N]g[+N] 

^ \ _ -r I ^ \ _ ^ ^ (4.14) 



if-/ ^[~N] \ ^- i ^[-N] ^[-2+N] 



lOn 



The claim in (i.ll) is a bit too strong and we will see further that it actually only holds inside some 



strips on the complex plane. 



10 



where S{e) := ]^5o(0-0j) 



(4.15) 



which gives for instance 



hm UN-a 
L—>oo 



n 




*Kn 



n 



N-a 



(^l 



'N/2+1] 



^[~N/2-l] 



N-a 



^l3N/2-a-l] ^l-N/2-a+l] 



"P 



[3Af/2-a-l] 



[-Af/2+a-l] 



n 



"p 



N-a 



*Kn 



gl+N/2 



^[_7V/2+a-l] ^[3N/2-a-l] n„ [^[-iV/2]] 

As a consequence, the large L hmit of equation ( |3.5D is 

-Lp. Ta,lTl% ^[-N/2-a+l] ^[-N/2-a+l] 



(4.16) 

(4.17) 
(4.18) 



YafliO) 



Ta+lflTa-lfl (^[-^/2+--l] ^\-N/2+a+l] 



Iia{[Sy- 



■N/2\-\ 



N/2] 



XCDD 



(4.19) 



where the factor Xcdd ^^^ added as another zero mode, necessary to transform the double 
poles and double zeroes of S^ into the simple ones ||2l|. We will also see in section 7.3.2 
that this factor arises in our Y-system formalism in a natural way. 

In particular, at a = 1, we get the ABA equation (periodicity condition): 



-iLsinh^ej ^ 



1 



,iL) 



(R) 



Q'j:;U{0j-i/2)Q'j^U9j-i/2) 



Xcnni0i)S{OjY Q^il^{6, + i/2) qS?1i(% + i/2) 



(4.20) 



which expresses the fact that Yi^q{9j + iN/A) + 1 = (here Ti^i was replaced by the single 
surviving, last term of ( [4.8D ). 

In conclusion, we have shown here that the Y system implies the familiar ABA equa- 



tions 1 25, 21 1 . In the next sections we will see how these ABA equations for the spectrum 
of PCF can be generalized to any finite size L. 



5. Expressions for the energy of excited states 



There exists no full proof procedure to generalize the formula ( |2.10| ) to the excited states. 
The analytic continuation of H with respect to the mass is difficult, if possible at all for a 
general state in our model. The procedure of [^, ^ claims that for the excited states a set of 
logarithmic poles (different for each state) appears under the integral in ( |2.10| ). From (|4.1l| ) 
and the ABA we know that at very large L, Tafi{9) ~ (f){9 — iN/A + ia/2), where (t){9) = 
WA9 — 9j) is a polynomial encoding all real roots^^ For finite L the roots 9j will be shifted 
and in general become complex. These exact Bethe roots, as opposed to the approximate 



^^We consider here for simplicity only the situation when the Bethe roots 9j are real in the asymptotic 
limit. The case when they occur in complex conjugated pairs should not be very different but at the moment 
we did not try do to it. 
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Figure 2: Analyticity of the integrand cosli(^6') log ((1 + yi,o)(l + ^2,0)) and manipulations with 
the contours when iV = 3 



ones given by ( 4.20| ), are defined by the exact Bethe equations Ta^i 



Ja) 



-iN/i + ia/2) =0. 

There is a whole family of such roots when a G [0,A^], because even though the two 
functions Tafl{9) and Ta+ifl{9 + i/2) have the same limit at large L, they do not necessarily 
have the same roots at finite size. Each of these roots also gives rise to two zeroes and two 
poles in the y- functions, namely, as we see from ( p.l5 ), l + Yafi{0j—iN/4L+ia/2^i/2) = 

and 1 + Ya±ifi{0"' — iN/4: + ia/2) = 00. Among these families of finite size Bethe roots, 

.(f) 



for even N , and 



.JV±1) 



for odd N. We 



we will actually restrict ourselves to the roots 6- ' .„.„.„ , 

will argue that only those ones will contribute as poles caught by an integration contour. 

Separating the logarithmic poles (where 1 + l^a cancels) in the contour integral ( 2.1C| ) 
should give a familiar contribution X^jCosh^^j to the energy of a finite L state. This 
appears to be the right, though not completely well understood and justified, answer for 
some models, including the PCF at A^ = 2 plj . 

For PCF at A^ > 2 this procedure encounters another difficulty: the zeros under the 
logarithm in ( |2.1[1| ) appear to correspond in general to complex Bethe roots 6j. We have 
to decide what is the right integration contour in ( 2.1[1| ) when this formula is applied to 
an excited state. We are not aware of any well justified procedure for fixing the contour 
but we shall try to guess it on the basis of our numerical observations and the symmetry 
considerations. 

In the rest of this section, we consider the formula for excited states of the U{1) sector 
- the one which corresponds to the wave function |^) having the the maximal value of total 
spin Sl = Sr = M/2 w.r.t. the SU{N)ji and SU{N)l symmetries. In this case J^ ' =0 
and there are no auxiliary roots in the ABA Q-functions ( |4.6| ) (all of them are equal to 1 
except (5]y = If, see ([4.6fj4.7| ). In what follows, we shall distinguish even and odd A^'s. 

5.1 Energy of state in the U{1) sector at odd A^'s 



It is believed that the energy of an excited state can be obtained from ( 2.10 ) by an analytic 
continuation in the parameter L. This continuation has the effect of appearance of new 
singularities of the integrand in the physical strip in ( p. 10 ) and a certain choice of the 
integration contour, enclosing some singularities of the integrand j^, Q|. How it happens 
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in each particular model or state is usually a rather complicated question. It implies the 
analysis of positions of these singularities at a finite L but the large L asymptotics often 
serves as an important guiding principle. 

Here we propose a formula for the energies of excited states in the U{1) sector which 
seems to work well for any odd N . It is based on our numerical and analytic observations, 
in particular for the A^ = 3 case. It reads as follows 

^(^) = -^ E S?Y T"'"'"'"' d^ cosh (1^) log (1 + Y.AO)) (5.1) 

SO that we have the straight integration contours parallel to the real axis and shifted by 

Let us explain the reason for such a choice of contours. First, let us note that in order 
to have a real energy from (^]^) we should impose the following property of Y-functions 
under complex-conjugation: Ya^s{(^) = lAr_a,s(^)- We will restrict ourselves to the gauges 
where this property is a consequence of the relation 



Ta,s{0) = TN-a,s{&) (5-2) 

For finite L, we will focus on the roots 6j defined^^ by Tn-i q{9j + i/4) = 0. 



rp+ rp- 



Due to the very definition of H-y„ n = ^r — °' ^°' — , each 6i gives rise to two zeros and 
poles. In particular, 1 + Yn-i q{9) has a zero and a pole^^ at respective positions 6j — i/4 
and 6j — i/A because r^_i {9j — i/4) = and Tiv-i q{9j — i/A) = 0. In the large L limit 

these zero and pole almost coincide since 6j is almost real. By complex conjugation, we 
can also say that 1 + Yn+i q{9) has a zero and a pole at respective positions 9j + i/4 and 
9j + i/4. 

This structure is illustrated for A^ = 3 in figure |2[ Prom the Liischer corrections^^, we 
can say that the pole occurs below the zero for 1 + Yi^ and vice versa for 1 + Y2fi, at least 
for roots with even momentum numbers^^. This is important to ensure the right answer if 
we want the contours to be straight. 

In (|5.l|) we chose the integration contour to pass, for the ^^^ and — |^'th term in the 
sum, between those zero and pole. Deforming the contour to the real axis and computing 



^^One of the advantages of this straight contour is that it can be easily implemented in numerics. We 
will see indeed that the Y functions can be most easily computed on exactly these lines. We will also see 
that the statement holds only for roots with even momentum number, but for odd momentum number, the 
(slightly modified) contour stays very close to this straight line. 

13 ( "-! ) 

In this section, we will denote 9j for 6- ^ , because the other types of finite size roots don't contribute. 
^^In addition to this zero and pole, 1 + Yn-i „(S) has another zero at 6j + 3i/4 and a pole at each root 
of T n-3 y, but this will not have any consequence in the contour argument. 



^^In section S.l, we detail how this is proved in the asymptotic limit. Our numerics suggests that it is 
still true at finite size, and even in the conformal limit. 

^®So that the contour will actually have to be slightly modified for roots having odd momentum number. 



This will be done in such a manner that (5.3) will stay true 
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Figure 3: Analyticity of the integrand for the Energy and choice of the contours for A^ = 4 



the contributions of the logarithmic poles enclosed by the contour during that deformation, 
one gets the following formula 



E{L) 



' N 



ti ^i^(F) 



+^E 



cos 



m- 



2N 



sm 



N 



d^ cosh ( —e ) log (1 + Ya,o{e)) 



sinh ( ^ (% - i/4)] - sinh (^ {9, + V4) 



(5.3) 



In the thermodynamic limit (L ^ 1), the Bethe roots 6j become real and the second 
line of ( ^.3| ) reduces to the asymptotic result ( p.7D , whereas the term in the first line appears 
to be 0(e-"^^). 

5.2 Energy of state in the U{1) sector at even A^'s 

When N is even, the corresponding contour cannot be chosen as a straight line. We will 
conjecture here the analogue of (|5.3|) to be simply 

V/ pai9)log{l + Ya,oie))d9 +y^coshiej) (5.4) 



E{L) 



m 

'n 



a=l --oo-^l+al 



where the roots 9j are defined by Tn_ q{9j) = 0, so that the second term in ( ^.4] ) is real due 
to the reality of Tn q. The corresponding contour is shown in figure!^ 

We should admit here that this formula for the masses at even A^ has a status of a 
natural conjecture. We have not enough of numerical, or analytic evidence to be 100% sure 
in it. It would be good to verify it at least for the mass gap at A^ = 4, numerically and by 
means of the Liischer corrections at large L. 



6. Hirota form of Y-system and wronskian solution 

For the principal chiral field, Ta^s is defined for a = 0,1, ... ,N, while 1^^^ is defined for 
a = 1,2, . . . ,N — 1. We can solve the Hirota finite difference equation ( p.l2| ) (and the 
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corresponding Y-system) with the appropriate boundary conditions using its integrabihty. 
Any solution of ( |2.12| ) is gauge equivalent to a solution where To^s{d) = Tq^^O — s|) and 
TN,siO) = T]\ifi{6 + s|). The most general solution can thus be expressed [^] as an A^ x A^ 
determinant, in terms of 2A unknown functions qj and qj ^'^: 

TaA^) = -(-l)"Det(c,,fc)i<,-fc<^ (6.1) 

where Cjk = q] \ -\ — (s + a + lH ^^)) ifA;<a 

(if N \\ 

and Cj^k = 1j\(^ + 7:\—s + a + l + — — 2fc|| ifA;>a 

At this point, qj is not necessarily the complex-conjugate of qj and the gauge freedom 
reduces to two independent functions g and g 

QjiO) ^ 9(0) ■ qm (6.2) 

WiG) ^ W) ■ rm (6.3) 

As an example of this determinant solution, the large L, spin chain limit solution 
corresponding to the U{1) sector state with roots 9i is easily identified, in the {R) (or (L)) 



gauge, by plugging the following definitions for qj into (6.1) 



q,{e) = q-{e) = ^— ^ (6.4) 

qN{0) = qWiO) = PooiO) (6.5) 

where {e^^o _ ^-^de^N-ip^ ^ ^ ^ J^^^g _ g^^ ^gg^ 

fc 

To see that this is the correct parameterization of the U{1) solution, first, we can convince 
ourselves that Ta-i = 0, Tq^q = V'(^ — z T^" ), and second, that it reproduces the Ti^g 



generated by (^-4.5) where all Qj(^)L<^ are set to 1. For the vacuum state Poo{d) 



Now we will explain how to generalize this large L solution to any finite L. 



7. Solution of the y-system for PCF at a finite size 

This section describes how to solve the y-system by reducing it to a finite number of 
non- linear integral equations (NLIEs), that can be solved in its turn by iterative numerical 
methods. 

We will focus on f7(l) sector states, although the method is in principle applicable to 
any excited state (see the discussion in subsection |9|) . 



^^The general solution doesn't assume that q's and g's are complex-conjugated. Nonetheless, our numerics 
has shown that at least for the states in the U{1) sector, it is sufficient to restrict ourselves to the solutions 
where q's and q's are complex conjugated. 
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7.1 Definition of tiie jump densities 

We propose here an ansatz for the finite size L solution by adding to the large L polynomial 
expressions (q^-S^) for g's certain finite L terms decreasing for 9 — > itcxD and exponentially 



small for L — > oo or — )• oo : The finite L g^-'s take thus the form 

q^{e) = + Fj{9) When j <N and Im{e) < (7.1) 

q-{9) = - + Fj{9) When j < N and Im{9) > (7.2) 

qj^{9) = P{9) + Fn{9) When Im{9) < (7.3) 

q^{9) = P{9) + F^{9) When Im{9) > (7.4) 



Where 



1 f°° fin) 
FJ9) = / ^^^dr? When/m(0)<O (7.5) 

2tTT J__^ 9 -T] 

I r°° fin) 

FA9) = / ^^^^dr? When/m(0)>O (7.6) 

2«vr J_^ 9 -T] 



and the polynomial P has the same degree as Poo = liinL-s>oo -P given by (|6.6f) ^^. As a 
consequence of these definitions, we have 

gjt+o! - ,[.-°] = hm g-[+^] - g]-^] = -/, (7.7) 

SO that fj is actually the jump between the functions Qj and qj on the real axis. 

Eqs.( |7.l[ - [7.4[ ) define qj only below the real axis and qJ above the real axis, so that the 
determinant only allows to compute Ta^s inside the strip lm(0) € [— '^^^ — ^^, — ^^" + 
^^]. We can already see that these strips are the minimal strips to compute the Y functions 



on the integration contour of equation ( p. 11 ) , and we will see that it enables to compute the 
exact^^ energy of states in the U{1) sector, at any length L. 

The jump densities fj{r]) are well defined on the real axis where they take only imagi- 
nary values, which follows from ( |7.7D and they are exponentially suppressed at large L or 
large cosh(-^ry), as can be inferred from ( |2.9|) . 

7.2 Relation to the analyticity of T functions 

From the ABA ( [4.19| ) and the finite size equation (3.5), we can see that 



N 
l + Yao > 1 When \Im{9)\ < — (7.8) 

' 6*— s>oo 4 

rp+ rp- 

which means that 1 + Y^ o = t^ — °'°t.°'° — has a proper behavior in this strip, being a 

' -t a + ljO-' a — 1,0 



meromorphic function regular at infinity. On the other hand, when \Im{9)\ = -r-, 1 + 1^^ 







^*The way we fix this polynomial will be explained in section 7.4, where the finite size Bethe equations 
are discussed. 
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up to the precision of our numerical procedure solving these NLIE's. 
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oscillates at Re{9) -^ oo, and it diverges when, e.g., \Im{9)\ e [^, ^]. By that reason we 



4 ' 4 



conclude that the analyticity strip of 1 + 1^0,0 = t °'°j-°;^ — is {6*, \lm{Q)\ < ^}. From this 



we can identify the strips where the asymptotics ( [4.11 ) hold for T^ 



^0,0 



> 991-^/2] ^i^gj^ j^(0) ^ ^ (7 9) 



Lcosh(^)->oo 4 



Ta,o\o<a<N T" "^ 9^'^"""^/'^ when \Imi9)\ < ^ + i (7.10) 

Lcosh(^)^oo 4 

These conditions ensure the proper analyticity of 1 + Ya^ = ^ — °'°y°'° — , and the boundaries 
of the analyticity strips of each 1 + Yafi are given by the boundaries of the analyticity strips 
of the corresponding T functions^*^. 

Now, since we know that the T functions are described by wronskian determinants, 
these analyticity strips suggest that 

Qj is analytic when lm(0) < 1/2 (7-12) 

qj is analytic when lm(0) > —1/2 (7.13) 

The fact that the analyticity strip for qj ends up at lui{6) = 1/2 is reflected for instance in 
the fact that T/v_i q is not analytic when lm(0) > N/A + 1/2, which explains that Y/v-i,o 
isn't analytic when lm(0) > N/i. 

The equations ( |7.12|j7T3| ) teach us that the analyticity domain is a bit bigger than 
what is necessary for (^-7^). It tells us that in the definitions (7.5,7^), the contour 



can be shifted up to ibi/2. In other words the functions fjiij) are analytic on the strip 
Im(r7) < 1/2. 

It is noteworthy that even with these contour deformations, the determinant expres- 



sions (fni-^) describe the function Ta,s{0) inside the strip Im((9) G]-^^^-f-l, -^^ + 
I + 1[, which is narrower than in equations ([7^, 7.11 ). But we will show that the relatively 



narrow strips given by this ansatz are sufficient to solve the y-system and compute the 
energies. 

7.3 Closed system of NLIEs 

The gauge freedom under the transformation (|6.2,6.3) can be used to impose Fi{d) = 



^i(^) = 0, which leaves only A^ — 1 independent densities to compute. A^ — 1 equations on 
this densities can be obtained by asking that the state is symmetric, i.e. that Ya^-s = ^,+s 
and inserting this requirement into the middle node y-system equation (p.lD^^, in the same 



manner as |11]. 



^"The analyticity strips for To.o and Tn,o can be chosen as half plane thanks to an appropriate gauge. 

^^As a consequence, we are solving the y-system under the two following constraints : Ya,s ~ 
g-ipa{ ) s,o X consta.s on the one hand, and T^ 2i ~ Ta,i one the other hand. This second constraint 
is specific to symmetric states (which includes the states in the U{1) sector), such that Ya,-3 = Ya,s- 
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Equation (|3.5D is then reduced to: 



2 I r^l+N-a]r^l-a]^ *^« 

On the other hand, we can write the determinant expression ( |6.1| ) for Ta~i (0 — i ^ ° ), 
and notice that Ca^j = g^" while Ca+ij = Qj ■ After subtracting one column from an- 
other in the determinant, there is a full column oiqj^^^^ —q, = —fj which is exponentially 
small and that explains the exponential suppression of Ta~i- Expanding the determinant 
w.r.t. these columns^'^ gives the following linear system relating /j's to T^ _is: 



N -2a 

fa,-l(y-i ^ 



j 

<iet{ck,i)k^j + det(cfc,/) k^j 
Where daj = {-ly ^^^-^ ^^^^ (7.16) 

These Ckj are the coefficients of the determinant ( |6.1| ) defining Ta^-i (9 — i 4^" ), and 



finally equations (7.14,7.15) can be recast into 




[+N/2\rj.[-N/2\ 
2^aajl,WJ/jl,WJ - e J„_i \ [2a~ZN/2]rj.[-N/2+2a\ 

This is a closed system of equations on {/j(^)}eGK because all coefficients daj, and all T's 
can be computed out of /j's through several convolutions. 

The solution of the y-system is therefore achieved by solving this system of A^ — 1 
equations on A — 1 densities. The simple inversion of the linear system ( [7.15 ) brings ( 7.17| ) 
into the form 

f,{d) = H^{{Mri)]k=2...N-i) (7.18) 

This Hj is a contraction mapping in some vicinity of fj = when L is sufficiently large since 
it leads to an exponentially small fj{0). This implies that in some vicinity of L = 00, the 
mapping Hj has a fixed point that can be numerically found through repeated iterations 
of i^. 

The way we solve y-system is therefore simply the iteration of (7.18) and a good news 



is that, at least for N = 3, even at very small L, this procedure seams to converge to a fix 



point of Hj , giving a complete solution of ( |7.18D and thus of the y-system 



The two terms in ( |7.16[ ) correspond to the fact that before expanding the determinant, we've added 

I subtracted columns to get one fuU column of q 
(which corresponds to a principal value in (7.5, 7. q)). 



^[+o]i [-0] 
and subtracted columns to get one full column of gj'^"' — q~ = —fj and the other one of — ^-^ — 
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7.3.1 Numerically workable form for the NLIE's 



One difficulty of tliis process is that in order to compute the right hand side of (7.17|) 



/ rp\+N/2]rj,[-N/2] \*Kn 

it is necessary to compute ( [ 2a-'3N/2 ] f-%/2+2a] ) • But we have already seen that the 

description it terms of g^-'s having cuts on the real axis allows to compute Tq^q{9) only when 
lm(0) < —-^ + i. So, for instance, To^q{9 + i^) cannot be computed when the spectral 
parameter 6 is real. The denominator can nonetheless be computed if in the convolutions 
we shift appropriately both the argument of the kernel and of the T-functions 



rj,[2a-3N/2]rp[-N/2+2a] 



i ] 



*Xiv'' 



rp[-N/2-N+a+l] ) I ^[+Ar/2+a-l] 



(7.19) 



since K]^{6) is regular when lm{9) G [— "^^y^, "^^2^] ^^^^ because a G [1, A^ — 1]. Eq.( |7.ig| ) 
simply reflects the fact that if k{6) is analytic for lm(0) E [0, b/2], then for real 6, [f^^') = 

y*{fcW)^ V/. 

The same idea, applied to the numerator, would give (To,o(0 + -^f)) ^ = {Tofi{9 — N^)) '^ . 
But the equality fails because Kn has a pole at i— g— • 
Instead, one can use the following relation 



.1 



*^^ To,o(0 -Ni + 



To,o{e + N-)] = "-^ ^ '-^ -^ (7.20) 

^0,0 ^0,0 ■■■-'0,0 

To,o(^ - ATI - 1) V O'O ; ^ ^ 

which simply uses the fact that for V/ regular, (Ilivi/])* ^ = /• This is true only up to a 
zero mode of Kn , discussed in the section 7.3.2| . 



Finally, the last factor in the eq.( [7.17D can be put into a numerically workable form by 
rewriting 

[+N/2]rj.[-N/2] \ *^^ 
rj,[2a-3N/2]rj,[-N/2+2a] 

r_Ar/2+l] [Ar/2-1] / rpl-3N/2-a+l]\*^N / rp[3N/2+N-a+l]\ *-^n 
^[-Af/2-1] rp[N/2+l] \ rj.[-3N/2+a-l] j I rp[3N/2-N+a-l] 



(7.22) 



This will help us to transforms the eq.( \l.l7 ), once the appropriate zero mode is added 



into a really closed system of NLIEs where the right hand side can indeed be computed by 
knowing the functions fj only on the real axis. 
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7-3.2 XcDD factor 

It is clear from the derivation of eq.( 7.17| ), as well as of the eqs.( [7.20 , 7.22| ) that they are 
fixed only up to a zero mode of the operator K^. A zero mode Z therefore has to be added 
to ( p7| ), to get 



^4,,(e)/,(0) = Ze 



-Lcosh{'^{B-i'i 



>=in(^) 



a~N/2] 



-N/2+1] rj.[N/2-l] 



'0,0 



T, 



Nfl 



rj,[-N/2-l] rj.[N/2 + l] 

-^0,0 ^Nfi 




^n 



p[3Ar/2+W-a+l] ' 

■Nfl 

p[ZN/2-N+a-l] 
'N,0 



[-N + a + l] 



(7.23) 



Such zero mode can include for instance the factors e '^°^ ( n > and Xcdd from equation 
( |4l9D . 

In the asymptotic limit (L — )■ oo), the zero modes in the equation (7.14) can be 



obtained by comparison with (4.19). In the same manner, the zero modes implicitly present 
in equation ( 7.22| ) can be computed in the asymptotic limit, by replacing Ta^s by their 



asymptotic values in terms of (p, see eq.(4.11), so that we can directly compute the zero 
mode Z when L — )• oo. We notice that, although both ( [7.17 ) and ( 7.22| ) are true up to a 
non-trivial zero mode, the zero mode in ( [7.23 ) happens to be^^ 



Z = l, 



(7.24) 



at least in the asymptotic limit. In the numerical solution of the y-system we therefore 
assume that Z = \ holds even at a finite size, i.e. that the analyticity structure of the 
zero modes is the same at finite L as at L — )• oo. We explicitly see that at L — )■ oo, Xcdd 
defined in ( [2.4[ ) is taken into account in ( [7.23 ). In addition we can check that at finite size, 
we obtain Y functions having simple poles only. 

7.4 Finite size Bethe equations 

Bethe equations emerge in this procedure as a regularity requirement on the jump densities 
/j's. Let us illustrate it for a general U{1) state in the SU{2>) case, and also show why 
these finite L analogues of Bethe equations are equivalent, at large L, to the ABA Bethe 
equations on the roots of (p. 

For such a state, the linear system ( 7.15| ) can be written as 




ri,_i(0-f/4) 



(7.25) 



^^ Actually the right hand side of ( 7.23) is itself defined up to a factor of e'" 
is defined up to a e^''^ ' ^ corresponding to the choice of the branch of the log. As a consequence, a more 



because any f^ = e^*'°s/ 



exact statement for (7.24) is Z — e^ n . The choice of k is done to reproduce (4.19), where the phase in 
.31) is chosen in such a manner that one particle at rest {9\ = 0) is a solution of the Bethe equation (4.20). 



For states with zero momentum (like the mass gap and the vacuum) , this extra phase can also be obtained 
by requiring a 6' — >■ —9 symmetry, which these states should exhibit. 
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where A = -^{q3 + qs) + Qs and S = 5(92 + 92) - ^2 • 

Inverting the matrix — — , some singularity could occur at the zeroes of the 
\-A-BJ 

function AB — AB, i.e. when the determinant is zero. If we want /j's to be regular, we 
need the numerator to vanish at the same 9 to cancel this pole. This gives the following 
finite size Bethe equation: 



One can notice that at such 9j the two conditions in the r.h.s. are equivalent. 

The claim that 9j are a finite size analogue of the Bethe roots is supported by the fact 
that at large L, the roots of AB—AB are precisely the Bethe roots. Indeed, at large L, B c^ i 
and A ~ P++ — P, giving AB — AB ~ i{P — P + P++ — P) = iip. Moreover, the second 



relation in the r.h.s. of ( [7.26 ) reduces then to the reality condition y' le-iU) ~ ~^- 



Using the leading-order large L expression of Yafi in terms of 5, eq. (f4.19|) , we get at large 
L 

Ti _i(^ - i/A) ^ ^-^^±|^e-^-«^(f (^-^/4)) (7.27) 

2(/90 

where 5(0) = \[Sl[e - 0,)Xez,z,(^ " 0,) (7.28) 

j 

Using the fact that (p{9i) = at all Bethe roots ^j, and dividing by the complex conjugate, 
the large L regularity requirement becomes 



{V f f 



ifS— {^++ys 



iLsinhi^ej) ^ _^ (7.29) 



Using the crossing relation, the left hand side becomes simply S{9j)e^ '^"^ ^^^i\ so that 
the finite size regularity condition stated above is equivalent at large L to the asymptotic 
Bethe equations ( [4.20| ). 



As a consequence, the iterative solution of the closed, finite size equations ( [7.23 ), should 



start from the expression (7.1,[f!4) where P = Poo is given in terms of the asymptotic Bethe 



roots by ( [6.6| ), and then at each iteration, this polynomial is updated in order to incorporate 
this regularity condition. 

A^ > 3 case : The same construction leads to finite size Bethe equations for any odd N . 
Like in equation (7.26) the number of regularity constraints at each zero of the determinant 



is apparently N — 1 but reduces to only one constraint: namely the cancellation of the 
/ ri,_i(0-i(A-2)/4) \ 

projection of \ to the kernel of the matrix djj defining the 

\TN-i,-i{9 + i{N-2)/^)j 
linear system ( [7.15| ). 
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0.001 0.01 0.1 1 10 L 

Figure 4: Energies of vacuum, of mass gap and of some excited states as functions of L, at A^ = 3 

This procedure for finite-size Bethe equations was described here for odd A^ and for 
states having real Bethe roots in the asymptotic hmit. The subtlety which arises when A^ 
is even, or for the states having, in the asymptotic limit, complex-conjugated complexes 
of Bethe roots, is that the zeroes of the determinant do not lie on the real axis but ap- 
proximately on M lb i/2. The above procedure can in principle be applied anyway, but its 
interpretation is less clear because the regularity condition is imposed at the very boundary 
of the analyticity strip. 

7.5 Numerical results 

As seen in the figure Q, this method allows to compute the energies of U{1) sector excited 
states for a wide range of lengths L. We can see that in the IR, L —^ oo limit, the energies 
of individual states behave mainly according to the number of "particles" forming the state 
- the number of the Bethe roots 6j: The vacuum energy tends to 0, while the energies of 
the states 8q, 9i and 02 tend to 1, and the energy of ^qo tends to 2). In the conformal limit 
-L — >■ we will see that the behavior is defined by the "particle's" mode numbers: The 
energy goes to -^( j^ — \- n) where n is the total momentum mode number. 

Numerical restrictions As the length L is decreasing, the algorithm looks worse and 
worse converging, and the densities become more and more peaked around the endpoints 
of the distribution. By choosing a small enough interpolation step (the densities fj are 
numerically defined by polynomial interpolation from a finite number of values), it was 
nevertheless possible to make the algorithm reasonably convergent for the considered states 
and lengths L, when N = 3. Decreasing further the interpolation step means increasing 
the computation time, and the necessary amount of memory, which puts a limit to our 
precision and the minimal length^^. 



^*We are aware of the fact that our precision for the energy levels is not sufficiently good. Our numerical 
results are still quite preliminary, we will try to improve them in the nearest future by using a better 
interpolation procedure and using more powerful computers. We definitely believe that our procedure can 
give a much better precision than these first results in the A'' = 3 case, and we hope also to tackle the A'^ > 4 
case. 
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L/sm(f) 


^vacuum 


Eeo 


Ee, 


Ee, 


Eea.Q 


60 


-(7) 10-24 


1.00000(0) 


1.00(7) 


1.28(6) 


2.00(3) 


20 


-5.1(6) 10-s 


1.00000(0) 


1.06(3) 


1.23(5) 


2.02(6) 


10 


-4.2(0) 10-4 


0.99(8) 


1.23(3) 


1.76(1) 


2.08(2) 


8 


-2.80(1) 10-3 


0.98(9) 


1.3(4) 


2.0(6) 


2.1(0) 


4 


-0.119(7) 


0.87(7) 


2.0(2) 


3.7(0) 


2.1(8) 


2 


-0.782(3) 


0.45(2) 


3.4(2) 


6.8(2) 


2.4(3) 


1 


-2.43(2) 


-0.59(2) 


6.2(0) 


13.(0) 


2.7(7) 


0.8 


-3.29(1) 


-1.15(7) 


7.5(8) 


16.(0) 


2.9(0) 


0.4 


-7.72(1) 


-4.21(0) 


9.(8) 


3.(1) 10^ 


3.2(7) 


0.2 


-16.9(0) 


-10.9(1) 


2.(8) IQi 


6.(1) 10^ 


3.2(6) 


0.1 


-35.(8) 


-25.(3) 


5.(7) 10^ 


1.(2) 102 


2.(0) 


0.04 


-94.(3) 


-71.(7) 


1.(4) 102 


3.(0) 102 


-5.(4) 


0.01 


-3.9(6) 10^ 


-3.2(2) 10^ 


5.(6) 10^ 


1.(1) 103 


-6.(1) 10^ 


0.004 


-1.0(1) 10^ 


-8.4(7) 10^ 


(1) 103 


(3) 103 


-(2) 102 


0.001 


-4.1(6) 10^ 


-3.5(8) 103 


(5) 103 


(1) 104 


-(1) 103 



Table 1: Numerical energies for several C/(l) sector states at A^ = 3 

Unfortunately, at A^ > 4 the calculations become heavier and with the size of inter- 
polation steps we can afford our algorithm becomes instable already for L of order ^ 1 
(which means we cannot really check the conformal limit for instance). At the moment we 
cannot say whether this has a physical meaning (like some symmetry breaking down, or 
some new type of singularity appearing) or whether it is just a numerical artifact, due to a 
poor numeric accuracy, or to the choice of the equations. For instance, it could be that the 
function we iterate over stops to be a contraction mapping but still has a fixed point, and 
maybe even that by rewriting slightly the functions it could become a contraction again, 
and extracting its fixed point would be possible by iterations. 



8. IR and UV limits 

8.1 Leading order results at large L 

The approach of this paper allows to compute the first exponential Liischer type correction, 
the so called Liischer correction, to the energy at large L, as we will show now on a few 
examples. 

8.1.1 Vacuum 

The large L behavior of vacuum is given by the condition that 



Yafl\ L^oo = (ra,l)2e-^P'" 
Liischer 



^.1) 



23 



where Ta^i is equal to the binomial coefficient . This is obtained from ( 4.19 ) where 



(f = 1, and can be plugged directly into ( p.lOf ) to get the energy to the leading order. By 
construction, the resulting energy fits well our numeric results when L is large enough^^. 

8.1.2 Mass Gap at A^ = 3 



When A^ = 3 it suffices to compute Yifi to get the energy, because l2,o = ^i,o- 
Moreover the previous analysis shows that 



Yio = e-^^^^'^^T 



70,0^2,0 5o(^-3ir v>{o+i) xcnnio-n: 



^-Lcosh(^.) (3^-5y 1___ l___ (33) 



-L cosh(4f e) y^^" ^4> 

+ if So{9-3lYXcnniO-3l 



that enables to compute at large L the leading order value of the integral term in (5.3). 

The analogous analysis for the second term in ( |5.3D is a bit more tricky, as it involves 
the position of the Bethe root. This position can be estimated by computing the densities 
to the leading order, to deduce the ffi'st correction to Ti^ in order to solve the equation 
ri,o(^o + i/4) = 0. 

For the mass gap, this root should be at the origin, up to exponential corrections 
in L. Moreover one can show^^ that Tifl{0 + i/4) ~ |/2(0) + 2/3(0) = ©(e--^^/^)^ ^j^jig 
r{ o(0+i/4) ~ i, so that Ti^o (6*0+2/4) = gives 6*0 ~ -|/2(0)-/3(0). Using the asymptotic 

expression for /j's (which can be extracted by keeping only the leading order in Ta^s and in 

je-^^/2r(-i)^r(2)2 
da j in the formula ( |7.17| ) ) , one gets ^o ~ ,l(i , so that the second term in 



( |5.3| ), which is sinh (^ (0o — V4)) — sinh (^ (^o + ^4)) can be computed at leading order 

That gives 

/32e-^^/2 3\ 

which is in very good agreement the numerical results, as can be seen in fig. 

Moreover, this expression (8.4) coincides exactly with the so-called /i-term |35, B3] 



which is known to dominate the finite-size corrections in the presence of the bound states. 

8.2 Conformal limit at L — > 

Let us start from the vacuum. At very small L, the effective coupling constant becomes 
very small e'^{L) ~ |,^\| and we can linearize the field on the group manifold in the vicinity 
of g{a, t) = I as g~^d^g ~ df^A, where A{a, r) is a Hermitian N x N traceless matrix field. 
The SU{N) PCF model should become a 2d CFT of A^^ — 1 massless bosons: R{L) is very 
big, the action will be 

^^The precision is over 10% when L > 4. 

^^These large L expressions are obtained by neglecting integral terms in the determinant expression of 
Ti.o. 
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^- = -Wr\ / ^*^^ 5^tr(a«A)2 + O {et{L)) . (8.5) 

In the ground state, the Casimir effect wih define the hmiting energy: Eq ~ — f^ + 
0(l/log'^(L^^)), with the central charge c = A^^ — 1, which gives -E'o^ — ""^^i^- 
The energies of excited states are 



fc=l o=l 

where n^ = (nj^ , n^ , • • • , n^ ) are the momentum numbers of particles constituting 

the state. We see that the small L asymptotics of our plots are well described by this 
formula: The vacuum, and the states Oq, do^, . . . have total momentum zero, and their 
energy satisfies ■^E{L) ~ — ""^ ^^ ' . On the other hand, a state like 9i has momentum 
number equal to 1 and ^E{L) ~ — ^^ ^~ ' + 1, etc. 

The qualitative behavior of the states ^Oi^ooi etc, at very small L's can be explained 
by the fact that the quantum fields are dominated by their zero modes. Since the momen- 
tum modes are not excited the field Y(a, r) does not depend on a. The action and the 
hamiltonian become: 

where the g(r) represents the coordinate of a material point (a top) on the group manifold, 
and J is the corresponding angular momentum on the group. The quantum mechanical 
spectrum of this system is well known: the quantum states are classified according to 
the irreducible representations of su{N) characterized by highest weight with components 
(mi > m2 >, • • • , > ?Tijv) usually represented by a Young tableaux A with N rows with the 
lengths m,j, j = 1, • • • ,N. The operator tr J-^ is nothing but the second Casimir operator 
with the well known eigenvalues, so that 

2^^^^-^'^^m^''^^^''-^ iviog(i/L) (^-^^ 

where r^ = m^ — 7vS7=i"^i- -^°^ instance, for a state with only M real roots in the 
asymptotic limit (and without self-conjugated complexes of roots), we have mi = M, 
mk>2 = 0, and 

^(F N-IM{M + Nl 

AI times 

This formula explains well the fact that the corresponding plots on fig.H| converge 
slowly, as inverse logarithm of L, to — (A^^ — 1)/12. 

The perturbative calculation of the mass gap [Eg^{L) — Eq{L)] for L <C 1 was done in 



1 37] and was compared with the numerical results following from the TBA approach in Q. 



We only recall that for the mass gap, in the logarithmic approximation. 
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E-L 
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Figure 5: Numeric energy of tire mass Gap (blue dots), compared to tlie analytic expression of 
E^^^°'P [green dotted line] and to the conformal limit Eq^^ ~ E^ac + I ^ ,___,_,,, , i }___„___, _,r^^ fred 
dashed line], where c — 12.3 is chosen to fit the data. 



9 L log(c/L) + i log(log(c/i)) 



L 

2^ 



[Ee,{L) - Eo{L)] 



N^-l 



1 



iV2 



+ k log log 



{L « 1) 



.10) 



logf ' 2 

which is in the perfect agreement with ( |8.8|) , as well as with our numerics. 

The formula (|8.8| ) also gives a prediction for state ^oo with two zero-momentum- 
particles, namely e''°(l)-e (l) ~ T' This result is consistent with our numerics when 
L is between approximately 10~^ and 1 (see the fig. |5|). For smaller L, although the 
algorithm seams to converge, the numeric precision is to poor because the integrand of 
(5J) is very peaked, and the compatibility with L — >■ predictions cannot be checked very 
accurately in this regime. 

Although our approach is based on adding some terms to the finite size solution, it 
reproduces correctly these conformal expressions, which shows that this description is not 
only accurate in some vicinity of L = oo, but even in the conformal limit where L is very 
small. This confirms that these terms where added in a sufficiently general manner to 
describe the relevant solutions of the y-system. 



9. Discussion 

We have presented here, on the example of the SU{N) principal chiral field model, a 
potentially powerful and rather general approach to the study of the finite volume spectrum 
of various integrable 1-1-1 dimensional sigma-models. The approach continues the ideas of 
|11] where the method was proposed on the example of the SU{2) PCF, but for N > 2 
the method has to be seriously reconsidered due to many new physical features w.r.t. the 
N = 2 case. In particular, the presence of the bound state particles and the non-reality of 
the Bethe roots at finite L show a few qualitatively new features within our approach. 

For virtually all integrable sigma models, the TBA-like approach of Al.Zamolodchikov 
can be summarized in a very universal system of functional equations, the y-system. The 
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Y-system is equivalent to the famous Hirota equation - the Master equation of integrability 
describing in this case the integrable discrete dynamics with respect to a pair of "represen- 
tational" variables, a, s and the spectral parameter (rapidity) 6. The boundary conditions 
for a, s are defined by the symmetry algebra of the model, whether as the analytic structure 
w.r.t. the 9 variable is in general the most complicated issue, really defining the dynamics 
of the model. It is plausible that in fact even the possible analyticity structures are greatly 
constrained by Hirota dynamics and by the symmetry algebra. It would be interesting 
to classify possible types of analyticity stemming from Hirota dynamics and some simple 
physical arguments (relativistic invariance, crossing, absence of certain singularities, etc.) 
related to the finite volume sigma models, similarly to the S-matrix bootstrap theory of 
Al.&A.Zamolodchikov valid only at infinite volume. This could lead to an interesting clas- 
sification of sigma models themselves and eventually to the discovery of new integrable 
models. 

In this paper, we managed to transform one such relativistic a-model, the SU{N) x 
SU{N) principal chiral field into a finite set of NLIEs, by solving its finite L Y-system 
in terms of wronskian determinants of a finite number of Q-functions and parameterizing 
these Q-functions by A^ — 1 densities correcting their large L asymptotics to any finite L. 



Generalizes the analytic and numerical results of |11| to A > 2, and we could numeri- 
cally check, at least when A^ = 2, 3, that this procedures solves the l^-system, and enables 
to compute energies for a wide range of lengths L, compatible with the UV conformal 
limit and the IR finite size (Liischer) corrections. On the way, we conjectured the natural 
generalization of the energy formula for excited states for the U{1) at finite L, to A > 2. 
This generalization appears to be unexpectedly non-trivial and looks different for even and 
odd N. The question of definition of the energy formula for excited states deserves a better 
understanding and hopefully the eventual derivation. 

The analysis was done for U{1) sector states, and it certainly can and should be 



generalized to any excited state, as was done in |11] for N = 2. To do this, one will 
have to understand the asymptotic terms, and the structure of zeroes. In particular some 
extra zeroes would appear, which might affect the way the energy is computed by contour 
manipulation. Apart from that, the main difference with U{1) sector would be that for 
non-symmetric states (i.e. when Ya^s 7^ ^,-s)) it would be necessary to introduce A — 1 
densities for the right wing and A^ — 1 densities for the left wing. One would have to write 

( |3.5| ) as 2(A — 1) different equations, by writing the left hand side either as \'fl) '^[r] or as 

'^2,0 ^0,0 

rp{L)rp[L) 

^(L) ^'(L) ■ But a similar approach based on the wronskian solution of Y-system should a 
priori still enable us to compute the energies of these states. 

An interesting problem which our approach might help to solve is the planar N —?■ oo 
limit in PCF. This model has a rich history of its comparison to QCD and it might provide 
an important example of the exactly solvable 2-|-l dimensional bosonic string theory, sim- 
ilarly to the matrix quantum mechanical model of the 1-|-1 dimensional, c = 1 non-critical 



string theory proposed and solved in |4C]. An exact and explicit solution for it was given 



in the case of infinite volume L but in the presence of specific "magnetic" fields |41] 
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As concerns the numerics, our algorithm converges very well for any length when 
A^ < 3, but for A > 4 it is very unstable for small enough length L, already at L < 1 
(which means for instance that we cannot really check the convergence to the conformal 
UV limit). Hopefully this instability has no direct physical meaning and is just a numerical 
artifact, due to a poor numeric accuracy, or to a bad choice of the iteration procedure for 
our NLIEs. 

We believe that this method of deriving a finite system of NLIEs for integrable sigma 
models is general and powerful enough to work for much more complicated cases of AdS/CFT 
correspondence, such as the superstring on the AdSs x S^ background dual to N=4 SYM 



theory, and the so called ABJM model where the y-system was already discovered [14, 38 



p9|| . The wronskian quasiclassical character |17] and even the full quantum solution of the 



Hirota dynamics for AdS5/CFT4 [^] are already available. It is left there to understand 
the very rich and complicated analyticity structure of Q-functions for short operators to 
complete the derivation of the AdS/CFT NLIE. 
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